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In the recent years, remarkable attention is drawn to Stirling engine due to noticeable advantages, for 
instance a lot of resources such as biomass, fossil fuels and solar energy can be applied as heat source. 
Great number of studies are conducted on Stirling engine and finite time thermo-economic is one of 
them. In the present study, the dimensionless thermo-economic objective function, thermal efficiency 
and dimensionless power output are optimized for a dish-Stirling system using finite time thermo-eco¬ 
nomic analysis and NSGA-II algorithm. Optimized answers are chosen from the results using three deci¬ 
sion-making methods. Error analysis is done to find out the error through investigation. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

At this moment it is very well known that the Finite Time Ther¬ 
modynamics has been criticized because did not take into account 
the internal irrevesibilities in the study of the real cycles [1-3]. Be¬ 
cause of that any optimization based on it would not get results 
very closed to the reality of Thermal Machines. In order to improve 
the situation Ibrahim et al. [4] introduced for the first time a “coef¬ 
ficient of irreversibility” which would take into account the inter¬ 
nal irreversibilities. This was done similar with the one introduced 
by Chen in [5] (in 2012) and implemented in their formulas but 
much earlier, in 1991. Unfortunately despite of the fact that this 
idea “improved the situation” of Thermodynamics with Finite 
Speed this could not get the Validation, drawing parallel with the 
computation with experimental results. Only the progresses made 
in the last 20 years in the field of Thermodynamics with Finite 
Speed [6—12] made it possible to get a Validation for 12 Stirling 
Machines [9] and 16 regimes of functioning. Based on the previous 
literature surveys [6-8] a very important Validation was obtained 
for Solar Stirling Engines [10-12]. 


* Corresponding authors. Addresses: Faculty of Mechanical Engineering, K.N. 
Toosi University, Tehran, Iran (M.H. Ahmadi), Institut de Recherche en Genie 
Chimique et Petrolier (IRGCP), Paris Cedex, France (A.H. Mohammadi). 

E-mail addresses: a.h.m@irgcp.fr, mohammadhosein.ahmadi@gmail.com 
(A.H. Mohammadi). 

0196-8904/$ - see front matter © 2013 Elsevier Ltd. All rights reserved. 
http://dx.doi.Org/10.1016/j.enconman.2013.05.031 


Performance analysis of heat engines applying finite-time ther¬ 
modynamics were begun by Chambadal [13] and Novikov [14] and 
were followed by Moran [3], In the first place, the performance of 
an endoreversible Carnot heat engine at maximum power output 
was scrutinized by Curzon and Ahlborn [15]. Many optimization 
researches on heat engines based on endoreversible and irrevers¬ 
ible models have been done with respect to finite-time and fi¬ 
nite-size constraints under various heat transfer modes, largely 
linear and non-linear ones, over the decades. Some optimization 
works could be found in literature survey represented by Bejan 
[16], Chen et al. [17] and Durmayaz et al. [18]. Typically, the opti¬ 
mization criteria are; the power output, thermal efficiency, entropy 
generation, the ecological benefit and thermo-economic objectives. 
The criteria were utilized by Sahin and Kodal [19], Kodal and Sahin 
[20] on the endoreversible and irreversible heat engines. Newly, 
the optimum operation conditions of an endoreversible solar-dri¬ 
ven heat engine with different heat transfer laws in the thermal 
couplings, operating at maximum ecological function conditions 
are studied by Barranco-Jim’enez et al. [21]. The thermoeconomics 
of an endoreversible heat engine model under maximum ecologi¬ 
cal conditions is investigated by Barranco-Jim’enez and Angulo- 
Brown [22], considering different heat transfer laws in a Novikov 
model. In thermoeconomic analysis of power plant models, an 
objective function is determined in terms of the cost involved in 
the performance of the power plant and a characteristic function 
such as power output [23] and ecological function. The thermoeco- 
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Nomenclature 



Ares 

absorber (recuperate) area 

f 

cyclic period, s 

^ app 

aperture area 

X 

temperature ratio of the Stirling engine 

c 

concentration ratio 

X 

vector of decision variables 

d 

deviation index 



di+ 

deviation or distance of ith solution from the ideal solu- 

Creek 

letter 


tion 

A 

ratio of volume during the regenerative processes (vol- 

di- 

deviation or distance of ith solution from the non-ideal 


umetric ratio) 


solution 

11 

thermal efficiency 

F 

dimensionless objective function 

8 

emissivity factor 

h 

l 

heat transfer coefficient, W K~ 4 or W m~ 2 K 1 
direct solar flux intensity, W irr 2 

s 

Stefan-Boltzman coefficient 

i 

j 

ith objective 
jth solution 

Subscripts 

n 

mole number of the working fluid, mol 

H 

absorber (heater) 

p 

dimensionless output power 

f 

convection at the high temperature side heat exchanger 

Q 

heat transfer, J 

L 

heat sink 

R 

universal gas constant, J moG 1 K 1 

k 

convection at the low temperature side heat exchanger 

S 

Entropy 

m 

entire solar dish Stirling system 

T 

temperature, K 

t 

Stirling engine 

W 

work, J 

0 

Ambient condition, optics 

V 

Volume 

4- f 

process states 


nomics of a Novikov power plant model is scrutinized by De Vos 
[24] in terms of the maximization of an objective function de¬ 
scribed as the quotient of the power output and the performing 
costs of the plant. 

Investigation of Solar Collector Design Parameters Effect onto 
Solar Stirling Engine Efficiency was studied by Ahmadi and Hos- 
seinzade [25], Intelligent approach to estimate solar Stirling heat 
engine power by implementation of evolutionary algorithm was 
extended by Ahmadi et al. [26]. Tlili investigated effects of regen¬ 
erating effectiveness and heat capacitance rate of external fluids 
in a heat source/ sink at maximum power and efficiency [27]. A 
mathematical model for the overall thermal efficiency of solar 
powered high temperature differential dish-Stirling engine with fi¬ 
nite heat transfer, irreversibility of regenerator, optimizing the ab¬ 
sorber temperature and corresponding thermal efficiency was 
evolved by Yaqi [28], 

Sharma et al. [29] studied the optimization of a solar-powered 
Stirling heat engine with the finite-time thermodynamics. 

The multi-objective optimization obstacle is a very tough pur¬ 
pose to be obtained which needs to satisfy a number of distinct 
and even contradictory objectives. At first, evolutionary algorithms 
(EA) were developed and applied in the 18th century to stochasti¬ 
cally solve problems of this generic class [30]. A logical alternative 
to a multi-objective problem is to interrogate a collection of solu¬ 
tions, which satisfies the objectives at a decent level without being 
overcomed by any other solutions [31]. Generally, Multi-objective 
optimization problems illustrate a large number of solutions so- 
called Pareto frontier, whose evaluated vectors represent the best 
trade-offs result. In this field, multi-objective optimization of dif¬ 
ferent themodynamic and energy systems have been inquired by 
researchers [32-39]. 

In the present investigation, multi-objective optimization is 
conducted to achieve the optimal dimensionless power and ther¬ 
mal efficiency of a dish-Stirling system with respect to thermo- 
economic parameter. Five decision variables such as temperature 
ratio, irreversibility factor, heat source temperature, hot working 
fluid temperature, ratio of heat transfer areas (heat sink to heat 
source) are considered. The presented analysis is great at designing 
and suitable design parameters can be obtained based on the anal¬ 
ysis. The paper is organized as follows: In Section 2, we present the 


thermodynamic and thermoeconomic description of the system, 
we apply the first and second laws of Thermodynamics to con¬ 
struct the objective functions, in Section 3 we describe the multi¬ 
objective optimization using evolutionary algorithms, in Section 4, 
we describe our decision-making methods implemented in our 
work, in Section 5 we present our main results and finally, in Sec¬ 
tion 6, we present our conclusions. 

2. Thermodynamic analysis of the system 

In solar-dish Stirling systems, mirrors of the parabolic shaped 
concentrator focuses the sun light on the focal point of the concen¬ 
trator where the hot end of the Stirling engine is located. Therefore, 



Fig. 1 . Schematic diagram of dish-Stirling system [28]. 
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the solar energy with a relatively high temperature is transferred 
to the hot side of heat exchanger of the Stirling engine. Fig. 1 illus¬ 
trates a schematic for a solar-dish Stirling engine connected to a 
solar dish concentrator. The solar-dish is equipped with a sun 
tracker which tracks the sun in order to have maximum solar en¬ 
ergy transfer to the engine when sun moves during the day. Hence, 
the solar energy is absorbed and transferred to the working fluid in 
the hot space of engine. 

As shown if Fig. 2, the Stirling cycle is consist of four processes. 
Process 1-2 is an isothermal process, in which the compressing 
working fluid at constant temperature, T c , rejects the heat to a heat 
sink at constant temperature, T L . Then the working fluid crosses 
over the regenerator and is warmed up to T h in an isochoric process 
2-3. In process 3-4, the working fluid expands at a constant tem¬ 
perature, T h , and receives the heat from the heat source at a con¬ 
stant temperature T H . Last process (4-1), is an isochoric cooling 
process, where the regenerator absorbs the heat from the working 
fluid. 

Actual useful heat gain of the dish collector, considering con¬ 
duction, convection and radiation losses is given by [26-29] as 
follows: 

flu = Muppflo - Arec [h(T H - T 0 ) + e5 (t 4 h - T 4 ) ] (1) 

where I is the direct solar flux intensity, A app is the collector aper¬ 
ture area, r/o = 0.85 is the collector optical efficiency, A rec is the ab¬ 
sorber area, h is conduction/convection heat transfer coefficient, 
T Ho „ is the (absorber temperature) heat source temperature, T 0 is 
the ambient temperature, e is emissivity factor of the collector 
and 5 is the Stefan’s constant. 

Thermal efficiency of the dish collector is obtained as [26- 
29]: 

fl s = ^ = flo - ^ [h(Tn - To) + es(t 4 h - To )] ( 2 ) 

where C is the collector concentrating ratio. 


Qp = [Mh(Th - T h )]t b (3) 

Q c = h c A L (T c - T L )t, (4) 


where q and t h are the times required for engine heat rejection and 
engine heat accumulation, respectively. T c and T h are the tempera¬ 
ture of the working fluid in the isothermal heat rejection and addi¬ 
tion process, respectively. A H denotes heat transfer area for heat 
exchanger of hot side and A L represents heat transfer area for the 
heat exchanger of the cold side in Striling engine. h h is high temper¬ 
ature side convection heat transfer coefficient and h c is low temper¬ 
ature side convection heat transfer coefficient. 

in the isothermal and endothermic process, the entropy term in 
the thermodynamic relation is: 

Tds - T ($), dv < 5 > 

The ideal gas equation is: 

PV = nRT (6) 


Substituting Eq. (6) into Eq. (5) yields: 

TdS = ^dV 

Integrating Eq. (7) over states 3 and 4 yields: 


w 


nRT 

IT 


dV 


Entropy is defined as: 
dS = 


(7) 

( 8 ) 

(9) 


Substituting Eq. (9) into Eq. (8) yields: 
Q 34 = Q„ = nRT h in (^j = nRT h in 2 
Let volumetric ratio 2 be defined as: 


( 10 ) 


2.1. The amounts of heat released by the heat source and absorbed by j = = Yl 

the heat sink V 3 V 2 


(ID 


The heat released between heat source and working fluid (Q/,), 
the heat absorbed between the working fluid and the heat sink 
(Qc) is obtained as follows [28]: 


In each cycle, variation in entropy in the working fluid inside a 
Stirling engine is obtained by: 



( 12 ) 


Solar 


Energy Convectio 



Fig. 2. Heat flows involved in a Stirling engine with solar collector. 


Since entropy is a state parameter it is independent of the inte¬ 
gration path, it should be noted that here, since processes are 
developed at the working fluid level ( T h , T c ), S gen is due only to 
internal irreversibility. Therefore, the net change in entropy be¬ 
tween the initial and final states of a complete cycle is identically 
zero: 


AC _Q* Q-c^c _ n 

AS — J T + Sgen — y T ®’ en — ^ 
From the Second Law of thermodynamics: 

Sgen 0 

The Clausius inequality is: 



(13) 

(14) 

(15) 


From Eqs. (13)-(15): 

(16) 

1 h 1 c 

Where the irreversibility parameter can be seen as a measure of 
the departure from the endoreversible case. Now, let [4,5]: 







M.H. Ahmadi et al./Energy Conversion and Management 73 (2013) 370-380 


373 


Oh 

T h 


CL_ 

<t> T e ’ 


<l> > 1 


(17) 


2.2. The conductive thermal bridging losses from heat source to the 
heat sink 


where </> is an irreversible property of a Stirling engine, such as ther¬ 
mal resistance, friction, heat loss. Consistent with the second law of 
thermodynamics, such irreversible factors, which cannot be ig¬ 
nored, bring about an entropy generation. Given an identical 
amount of heat transfer Qh, for an endoreversible heat engine <f>) 1, 
Qc is expressed as: 


Relating Eqs. (17) and (18) yields: 

Qc = <t>Q?c V 


(18) 

(19) 


In a Stirling engine, the heat that is released to low-temperature 
thermal storage (cold chamber) in a reversible isothermal process 
yield: 


Q_ r r = nRT c InQQ (20) 

Back-substituting Eq. (20) into Eq. (19) yields: 

Q c = <fmRT c In = cj>nRT c In 7 (21) 

Equating Eq. (3) with Eq. (10) yields: 


[hhA H (T H - T h )]t h = nRT/, In A 

Time of heat transferred process 3-4 is given by: 

f nRT h \nA 

th = Mh(T h - T h ) 

Equating Eq. (4) with Eq. (21) yields: 

h c A L (T c - T L )ti = (j>nRT c In A 

Time of heat transferred process 1-2 is given by: 


( 22 ) 

(23) 


(24) 


< pnRT c In A 
hcA L (T c - T l ) 


(25) 


Irreversibility of the finite-rate heat transfer causes the follow¬ 
ing fact that the time of the regenerative processes is not insignif¬ 
icant in comparison to the two isothermal processes [28]. To 
compute the time of the regenerative processes, the temperature 
of the working fluid in the regenerative processes is deemed as a 
function of time which is discerned by [5,26-29]: 


This value proportional to the temperature difference of the 
heat source and heat sink and the cycle time is obtained as follows 
[28]: 

Qo = fco {Th - T L )t (30) 

The net heat released from the heat source (Qh) and the net heat 
absorbed by the heat sink (Qi) are obtained as follows: 

Qh = Qti + Qo (31) 

Qi = Qc + Qo (32) 


2.3. Power and thermal efficiency 


Considering the cyclic period of the Stirling engine, the output 
power, the thermal efficiency of the engine are given by: 


p , = W = Q H — Q_l 


It 


t t 

Qh-Ql 
Qh 


Substituting Eqs. (22)—(32) into Eqs. (33) and (34) we have, 


P' = 


Th - <t>T c 


It 


h„A„(T H -T h ) + A R A H hJ(T c —T L ) + F 1 (Th Tcj) 

_ T h ~4>T _ 

T h + k 0 (T H - T l ) [ Mh( ^_ T|i) + ArA J*t c _ Tl) + F , (T h - T c ) 


(33) 

(34) 

(35) 

(36) 


Where = SRTiTt ^ 

For the sake of convenience, a new parameter x = ^ is intro¬ 
duced into Eqs. (35) and (36), then we have: 


P' = 


ft 


h h A H {T„-T h 


(1 ~ 4>x) 

+ A R A H h*xT„-T L ) + (1 -*)) 

(1 - tfx) 


Th + k 0 (T H 7l)[ Mh( ’ h 7- h) + ARAH/tcxTh-TiJ -I- f 1 (1 X)] 


(37) 

(38) 


where the A R represents heat transfer area ratio . 

The thermal efficiency of the dish-Stirling system is product of 
the thermal efficiency of the collector and the thermal efficiency of 
the Stirling engine [28,29]: 

rim = >1t1s (39) 


^ = ±JW, (26) 

where M is the proportionality constant which is independent of 
the temperature difference and depends only on the property of 
the regenerative material, called regenerative time constant and 
the ±sign belong to the heating (i= 1) and cooling (i = 2) processes 
respectively [5,26-29]. 


2.4. Thermo-economic model 

In the present inquiry, the objective function is determined as 
the power output per unit investment cost since there is no fuel 
consumption cost in a solar-driven heat engine [19,20,23]. For 
optimizing power output per unit total cost, the objective function 
is revealed as follows: 


J M, 

, _Th-T c 
U M 2 

The cyclic period 

Using Eqs. (22)—(28), we get that the cyclic period t is: 


(27) 

(28) 


t = 


nRT/, In A 


4>nRT c In A 


h h A H (T H - T h ) h c A L (T c - T L ) 


+ 


1 1 

m + M 2 


(Th ~ T c ) 


(29) 


F’ = r (40) 

where C ai refers to the annual investment cost. The investment cost 
of the plant is assumed to be proportional to the size of the plant. 
The size of the plant can be taken proportional to the total heat 
transfer area. Thus, the annual investment cost of the system can 
be given as [18-23] 

Cai = (oAh + Ml) 


(41) 
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Where the investment cost proportionality coefficients for hot- 
and cold-sides are (a) and (b), respectively, equal to the capital 
recovery factor times investment cost per unit heat transfer area, 
and the dimension is ncu/(yr m 2 ). 

Substituting Eqs. (37) and (41) into Eq. (40), we obtain 


F' = 


(1 - <t>x) 


1 -/ 


1 + 


h(TH~T h ) 


fjlX ^1 t 

' A R h c (xT h -T L ) 


+ W l-x)(l + 



(42) 


where / is relative investment cost parameter of the hot-side heat 
exchanger and defined as [18-23] 


/ 


a 

a + b 


(43) 


By using Eqs. (37), (39) and (42) one can obtain the dimension¬ 
less thermo-economic objective function F = ^, dimensionless 
power output P = hAuTL and thermal efficiency of the dish-Stirling 
system ij m , respectively, as 


F = 


(1 - <t>x) 


1+77 Mb 1+ 74 


p = 


Tl r (Th W + T h -T L ) ' + FjflhAh ^ x ) + (Y)^®) 

(1 - <t>x) 


Fl (ph-W + A R h c tX-T L ) + Fi ^AhCI *)) 


(44) 

(45) 


Pareto optimal if there is no other possible solution Y that would 
decrease some objective functions without concluding a simulta¬ 
neous increase in at least one other objective function. The collec¬ 
tion of all plausible non-dominated solutions in X is leaded to the 
“Pareto optimal set”, and for a certain Pareto optimal collection, 
the corresponding objective function values in the objective space 
are denominated the “Pareto optimal frontier” [31], 

Fig. 3 illustrates the objective functions space for an optimiza¬ 
tion problem with two objective functions and f 2 . The feasible 
are is insight and infeasible area is outsight of the curve. Each point 
in the feasible area denotes a solution. It is obvious that the values 
of functions/! and/ 2 . for point M are lower than the corresponding 
values of point j. Therefore point M dominates point J or point M is 
better than point J. In the same way points L and N dominate M. But 
points / and I< neither overcome M nor are overcome by M. Thus, 
only thee points that are placed in the left-down parts of M would 
overcome M. There is not any point in the feasible space placed in 
the left-down part of point R. Thus R is not overcome by the other 
points in the feasible space and hence R is a Pareto optimal. This is 
true for points P, A, R, E, T, 0 and all of the other points placed at the 
bold curve denoted as “The Pareto Optimal Frontier” in Fig. 3. In 
the feasible space, the minimum value of is for P and the mini¬ 
mum value of f 2 is for point 0. Thus P and 0 are the solutions for 
one objective optimization issues whose objective functions are 
/i and f 2 , respectively. The other points on the Pareto frontier are 
also optimal, but a decision-making procedure is needed to find 
which of them should be selected. 

3.2. Optimization via EA 


>1m = {'/o - ^ [h(T H - To) + 8 S{T 4 h - Tq)] j 

x _ (1 - fro _ 

Td + MFh - T l ) [ }lh ^ + ARA J^_ TL) + Fi(l - x)] 

(46) 

3. Multi-objective optimization with evolutionary algorithms 

3.1. General concepts of the multi-objective optimization 

Imagine a decision-maker that would like to optimize objec¬ 
tives which are non-commensurable. The decision-maker does 
not know the periority of the objectives. Moreover, all objectives 
are of the minimization type. For a minimization kind of objectives, 
it can be changed to the maximization type by multiplying nega¬ 
tive one. The procedure of a minimization multi-objective decision 
problem with objectives is represented as follows: 

Offer an non-dimensional decision variable vector x = {xi,... ,x n } 
in the solution area X. Detect a vector x* which minimizes a collec¬ 
tion of K objective functions /(x) = {/i(x*),.. .,/k(x*)}. The solution 
space X is limited by a number of constraints, such as gj(x*) = bj 
for j= 1 ,...,m, and bounds on the decision variables [31], Totally, 
it is not possible to minimize all the K objective functions at once 
with one solution vector X which means in real-world issues, 
objectives at issues disagree with each other. Thus, optimizing X 
with concerning a uni-objective mostly propels to unreasonable 
outcome with respect to the other objectives. Thus a novel impli¬ 
cation, so-called the “Pareto optimum solution”, is employed in 
multi-objective optimization quandaries. A practical solution X is 
denominated “Pareto optimal" if there is no other plausible solu¬ 
tion Y that prevails solution X. Therefore, a feasible solution Y is 
entitled to overcome another conceivable solution X, if and only 
if,/)(V) 7/;(X) for i = 1,.. ,,k and fj(Y) < /(X) for at least one objective 
function which denotes that a possible vector X is denominated 


In the present investigation, the Pareto frontier is achieved 
using the Genetic Algorithm (GA) which is identified as a field of 
evolutionary algorithm. Genetic Algorithms were evolved by John 
Holland in the 1960s to import the mechanisms of natural adapta¬ 
tion into computer algorithms and numerical optimization [33], 
They are applied as a computer simulation in a way that a popula¬ 
tion of abstract representations (called chromosomes or the geno¬ 
type of the genome) of candidate solutions (called individuals, 
creatures, or pheno-types) to an optimization problem develops 
toward better solutions. It commences from a population of ran¬ 
domly produced exclusives and occurs in generations. In every 
generation, the fitness of each individual is measured; multiple 
individuals are stochastically opted from the current population, 
and altered to form a new population. The new group is employed 



Fig. 3. Schematic of the objectives space. 
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in the next iteration. Generally, the algorithm finishes when either 
a maximum number of generations have been produced, or a sat¬ 
isfactory fitness level has been obtained for the population. If the 
algorithm has ended owing to a maximum number of generations, 
a satisfactory solution might have been reached or might not have 
been accomplished. In Genetic Algorithms, a candidate solution to 
a problem is named a chromosome, and the evolutionary viability 
of each chromosome is denoted by a fitness function. This method 
can be applied to optimize non-linear problems [31,34). 

Multi-objective evolutionary algorithms (MOEAs) have been 
evolved in the recent decades by a lot of tests on complex mathe¬ 
matical problems. On real-world it has shown that they can omit 
the obstacles of classical methods [31,34]. The structure of the 
MOEA exploited in this investigation is depicted in Fig. 4 [33]. 
The real values of decision variables are applied instead of their 
binary coded. 

3.3. Objective functions, decision variables and constraints 

Three objective functions for this study are the dimensionless 
thermo-economic objective function (F), the dimensionless power 
output (P) and the thermal efficiency of dish-Stirling system (i/ m ) 
denoted by Eqs. (44)-(46), respectively. 

In this paper five decision variables have been considered as 
follow: 

tj> Internal irreversibility parameter 

Ar Heat transfer area ratio 

x Temperature ratio 

Tn Temperature of heat source 

T h Temperature of the working fluid in the high temperature 
isothermal process 


The objective functions with respect to following constraints 


have been solved: 


0.45 < x < 0.7 

(47) 

1100 1400 K 

(48) 

0 > 1 

(49) 

850 < T h < 1000 K 

(50) 

0.25 <A r < 10 

(51) 


4. Decision-making in the multi-objective optimization 


optimization. At first, dimensions and scales of objectives space 
must be unified due to the possibility of different dimensions of 
objectives in multi-objective optimization issues. In turn, objec¬ 
tives vectors must be non-dimensionalized before the decision¬ 
making procedure. Some methods of non-dimensionalization 
including linear, Euclidian, and fuzzy can be employed in decision 
making. 

• Linear non-dimensionalization 


Consider the matrix of objectives at various points of the Pareto 
frontier is denoted by Fy where i is index for each point on the Par¬ 
eto frontier and j is the index for each objective in objectives space. 
Therefore a non-dimensionalized objective, FJJ-, is defined as, 


F 

F" =-L— For maximizing objective 

1 max(F,j) 

F 

F". =-For minimalizing objective 

max(l/ r y) 


(52a) 

(52b) 


• Euclidian non-dimensionalization 


In this method, a non-dimensionalized objective, FjJ, is defined as, 

F 

F". = . lJ . = For minimizing and maximizing objectives 

(53) 

• Fuzzy non-dimensionalization 
In this method, a non-dimensionalized objective, Fjj, is defined 


as, 




F n = 

V 

Fy - min(F,j) 

For maximizing objectives 

(54a) 

max(F,j) - min(F,j) 

F n = 

V 

max(F 0 ') - Fy 

For minimizing objectives 

(54b) 

max(F,j) - min(F,j) 


In the present investigation, some well-known and common 
kinds of decision-making procedures including the fuzzy Bell- 
man-Zadeh, L1NMAP and TOPSIS methods are applied in parallel 
and the final optimal solution has been opted based on real-world 
engineering experience and criteria among solutions which are 
recommended by the mentioned ways. The Bellman-Zadeh proce¬ 
dure implements the fuzzy non-dimensionalization while the 
other ways (L1NMAP and TOPSIS) use Euclidian non-dimensional¬ 
ization. The algorithms are introduced in the following. 


In multi-objective issues, a procedure of decision-making to opt 
the optimal solution from existing results is needed. There are lots 
of procedures for decision-making which can be applied to pick the 
final solution from the Pareto frontier that is achieved through 


4.J. Bellman-Zadeh decision-making method 

When using the Bellman-Zadeh approach, each FfX) is replaced 
by a fuzzy objective function or a fuzzy set, 



Fig. 4. Scheme for the multi-objective evolutionary algorithm used in the present study. 
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Aj = {X,h a .(X)} X e L. j= 1,2,... ,k (55) 

where n A .(X) is a membership function of Aj [40]. 

A final decision is defined by the Bellman and Zadeh model as 
the intersection of all fuzzy criteria and constraints and is repre¬ 
sented by its membership function. A fuzzy solution D with setting 
up the fuzzy sets (Eq. (56)) is turned out as a result of the intersec¬ 
tion D = rij_ | Aj with a membership function, 

/< D (X) = n'L 1/ u A (X) = min/(.(x), XeL (56) 

1 J = 1 

Using Eq. (54a), it is possible to obtain the solution proving the 
maximum degree as follows: 

max p D (X) = max min n A . (x) (57) 

X° = arg max min f.i A (x) (58) 

XeL j=l-.k > 


Ha( x ) 


1 - 


GjW-G^ 


if C,(x) > G™ ax 

if G™ x < Gj(X) < G™ x + dj 

if G,(x) < G” ax 


( 63 ) 


A final decision is determined as the intersection of all fuzzy cri¬ 
teria and constraints represented by its membership functions. 
This problem is reduced to the standard nonlinear programming 
problems: to find such values of X and k that maximize k subject to 


X < fi F ., i = 1,2,..., n 
Ucj, j = 1)2,... ,k 


The solution of the multi-criteria problem discloses the mean¬ 
ing of the optimality operator and depends on the decision-makers 
experience and problem understanding. 


4.2. LINMAP decision-making method 


To obtain Eq. (57), it is necessary to build membership functions 

H A ,xy.j = 1_, k, reflecting a degree of achieving "own” optimal by 

the corresponding F/X), X s I; j = 1. k. This is satisfied by the use 

of the membership functions [41], The membership function of 
objectives and constraints, linear or nonlinear, can be chosen 
depending on the context of problem. One of possible fuzzy convo¬ 
lution schemes is presented below [42]. 

Initial approximation for X-vector is chosen. Maximum (mini¬ 
mum) values for each criterion F/X) are established via scalar max¬ 
imization (minimization). Results are denoted as “ideal” points 

{x°J = l....,m}. 

The matrix table {Tj, where the diagonal elements are “ideal” 
points, is defined as follows: 



Maximum and minimum bounds for the criteria are defined: 


F™ n = min F,(X °), i = 1. n (60a) 

j J 

F™ x = max Fj(Xj), i = 1,..., n (60b) 


The membership functions are assumed for all fuzzy goals as 
follows. 

For minimized objective functions: 


MX) 


0 

pmaxp, 
pmax pmin 

1 


if Fj(x) > F™ x 
if F™ n < Fj s; F” 
if F f (x) < F™" 


For maximized objective functions: 

1 if Fj(x) > F™ ax 

F, ~ fr if Ff" < F, £ F™* 


MX) 


pmax _pmin 

0 


if Fj(x) < F™ n 
Fuzzy constraints are formulated: 

Gj(X) < G J max + d h j = l,2,...,k 


(61a) 


(61b) 


(62) 


where dj is a subjective parameter that denotes a distance of 
admissible displacement for the bound Gj max of the jth constraint. 
Corresponding membership functions are defined in following 
manner: 


The point that every objective is optimized irrespective to satis¬ 
faction of the others is an ideal point on the Pareto frontier. It is 
obvious that having every objective in the optimal condition is 
not possible for the multi-objective optimization and it can achieve 
in a single-objective optimization problem. In the result, the ideal 
point is not placed on the Pareto frontier. In the LINMAP procedure, 
after Euclidian non-dimensionalization of all existing objectives, 
the distance of every solution on the Pareto frontier from the ideal 
point marked by d,+ is defined as: 

di + = ( 65 ) 

where denotes the number of objective while stand for each solu¬ 
tion on the Pareto frontier (i= l,2,...,m). In Eq. (65), F , f eal is the 
ideal value for jth objective obtained in a single-objective optimiza¬ 
tion. In LINMAP method, the solution with minimum distance from 
ideal point is selected as a final desired optimal solution, hence, i in¬ 
dex for a final solution, ip nal is, 

ifimi = i 6 min(d, + ); t=f,2,...,m (66) 


4.3. TOPSIS decision-making method 


In this method beside the ideal point a non-ideal point is de¬ 
fined. The non-ideal point is the ordinate in objectives space in 
which each objective has its worst value. Therefore, beside the 
solution distance from ideal point, d,+, the solution distance from 
non-ideal point denoted by d,_ is implemented as a criterion for 
selection of the final solution. Hence, 



(67) 


In continuing the TOPSIS method a parameter is defined as follows, 


Cl,= 


d, 

d i+ + dj_ 


( 68 ) 


In TOPSIS method a solution with minimum C/, is selected as a de¬ 
sired final solution, therefore, if ip na i is an index of the final selected 
solution, we have, 

ijinai = i e max(C/,); i= 1,2,...,m (69) 


5. Results and discussion 

The dimensionless thermo-economic objective function (F), 
dimensionless output power (P) and thermal efficiency of the so¬ 
lar-dish Stirling System (f/ m ) are maximized simultaneously using 
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Fig. 5. Pareto optimal frontier in objectives’ space. 



Fig. 6. Pareto optimal frontier in objectives' space (iy m - P). 


the multi-objective optimizing method which works based on the 
NSGA-II algorithm, in this regard, optimization is performed with 
objective functions that are expressed by Eqs. (44)—(46) and con¬ 
straints which are expressed with Eqs. (47)—(51 ). Design parame¬ 
ters (decision variables) of optimization are the temperature of 
the working fluid in the high temperature isothermal process 
(T h ), temperature ratio (x), heat transfer area ratio {A R ), internal 
irreversibility parameter (<■/>) and the absorber temperature (T H) . 

in order to have consistency with previous works, specifications 
of the solar-dish Stirling engine are considered as follows [28]: 



Fig. 7. Pareto optimal frontier in objectives’ space (P - F). 

h f = h k = 200 W K 1 itt 2 , / = 0.7, C = 1300, 

<5 = 5.67 x 10~ 8 W nr 2 IC 4 , T L = 320 K, T 0 = 300 K, 
h = 20 Wirr 2 K _1 , J= 1000 Witt 2 , 

(Hl + S5)- 2xlrSsrl - s - 43J ™ rl 
n = 1 mol, 7 = 2, e = 0.9, k 0 = 2.5 WIC 1 , i/ 0 = 0.85 

Fig. 5 shows optimal results for P, ij m and F objective functions. 
Optimal results based on three decision making methods are 
shown and ideal and non-ideal consequences are indicated. 
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Fig. 8. Pareto optimal frontier in objectives' space (F - i; m ). 

Optimal solutions population in 0.45 ^P^ 0.55, 0.38 ^ r\ m 
■C: 0.41 and 0.15 19 can be seen in Fig. 5. Therefore, optimal 


answers are probably in this region and as it is mentioned, the re¬ 
sults of Fuzzy Bellman-zade, Linmap and Topsis decision-making 
methods are in this region. Ideal solution is in coordination 
(0.42,0.55,0.205) which are F, P and ij m , respectively. Non-ideal 
solution is in coordination (0.365,0.345,0.08). The chosen results 
through Topsis and Linmap methods are the closest to ideal solu¬ 
tion. With respect to 0.35-40% thermal efficiency of dish-Stirling 
system, optimal solutions region is appropriate. 

Optimal results for P and //„, objective functions are illustrated 
in Fig. 6 and obtained temperatures through decision making 
methods are displayed. 

Optimal solution population for two objective optimization 
(P - t]m) are in 0.392 ^ P C 0.524 and 0.397 0.421. The ideal 

and non-ideal solutions are in coordination (0.421,0.524) and 
(0.397,0.392), respectively. 

Fig. 7 demonstrates optimal results for P and F objective func¬ 
tions which are in 037 ^ P C 0.545 and 0.122 ^ F ^ 0.205. Non-lin¬ 
ear relationship between P and F can be seen. Ideal and non-ideal 
solution are (0.545,0.205) and (0.37,0.122), respectively. 

Optimal results for ij m and F objective functions are represented 
in Fig. 8. 

The results of three-objective optimization based on Fuzzy, LIN¬ 
MAP and TOPSIS decision making methods are shown in Table la. 
In this table, results for F, P and ;? m objective functions based on 
five decision variables are categorized. Temperature ratio (x) is ob¬ 
tained between 0.45 and 0.46 which is in acceptable range due to 
previous studies [16,17]. The achieved values for < j> are close to 1. 
4> = \ represents Carnot reversible engine. Table la demonstrates 


Table la 

Decision making of multi-objective optimal solutions. 


Optimal solution 

Decision variables 





Objectives 



X 

Th 

</> 

T h 

Ar 

Y]m 

P 

F 

TOPSIS 

0.45466 

1371.415 

1.002271 

920.6022 

3.175778 

0.384027 

0.436108 

0.184709 

LINMAP 

0.459212 

1371.466 

1.001516 

919.8346 

3.671255 

0.381995 

0.454229 

0.176510 

Fuzzy 

0.451457 

1329.38 

1.000984 

913.8662 

4.34813 

0.396756 

0.449328 

0.156916 


Table lb 

Decision making of multi-objective optimal solutions (ij m - P). 


Optimal solution 

Decision variables 




Objectives 


X 

Th 


Th 

Ar 

rj m 

P 

TOPSIS 

0.45004 

1340.77 

1.000153 

896.0217 

8.012228 

0.398154 

0.519160 

LINMAP 

0.450119 

1335.903 

1.000144 

896.0208 

8.016084 

0.399064 

0.515551 

Fuzzy 

0.450056 

1273.835 

1.000111 

896.0406 

8.005058 

0.410305 

0.466470 


Table lc 

Decision making of multi-objective optimal solutions (P - F). 


Optimal solution 

Decision variables 




Objectives 


X 

Th 

<i> 

Th 

Ar 

P 

F 

TOPSIS 

0.465289 

1382.388 

1.000085 

912.5242 

3.169516 

0.442871 

0.187788 

LINMAP 

0.461091 

1382.428 

1.000114 

911.3326 

3.555935 

0.457653 

0.181323 

Fuzzy 

0.455245 

1382.369 

1.000028 

911.1519 

4.072863 

0.474479 

0.172820 


Table Id 

Decision making of multi-objective optimal solutions p/ m - F). 


Optimal solution 

Decision variables 





Objective functions 


X 

Th 


Th 

Ar 

tJm 

F 

TOPSIS 

0.450132 

1373.37 

1.000107 

948.427 

1.798122 

0.383051 

0.204405 

LINMAP 

0.450049 

1357.56 

1.000065 

947.5719 

1.807338 

0.386469 

0.200713 

Fuzzy 

0.450135 

1262.96 

1.000066 

918.7098 

1.835071 

0.403693 

0.176783 
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Table 2 

Error analysis based on the mean absolute percent error (MAPE) method. 


Decision making method 

LINMAP 


TOPSIS 


Fuzzy 

Objectives 

Win 

p 

F 

y]m 

p 

F 

>lm P F 

Max error % 

1 

5 

6 

1 

4 

2 

1 10 10 

Average error % 

0 

1 

3 

1 

2 

1 

1 5 5 


that the efficiency of dish-Stirling system is between 38.4% and 
39.67% which is acceptable because of previous works [28,29]. 

The results of two-objective optimization (P and ;; m ) based on 
three mentioned decision making methods are represented in 
Table lb. The consequences of two-objective optimization for 
(F P) and objective functions are illustrated in Tables lc 

and Id, respectively. 

Mean absolute percentage error was carried out in error analy¬ 
sis of implemented methods. To assess this goal, 30 runs of each 
approach were performed to get final solution by bellman-zadeh, 
linmap and Topsis decision making methods. First row of Table 2 
demonstrate maximum absolute percentage error (MAAE) of three 
decision making methods. Moreover, second row of Table 2 report 
mean absolute percentage error (MAPE) of mentioned methods. Fi¬ 
nally it worth bearing in mind that, all of the above procedure en¬ 
coded in the MATLAB® software. 

6. Conclusions 

Throughout current research work enormous efforts have been 
made to figure out optimum condition of the three addressed 
objective function such as dimensionless output power, thermal 
efficiency of dish-Stirling system and dimensionless thermo-eco¬ 
nomic parameter while through developed multi-objective optimi¬ 
zation approach five different variables like as heat source 
temperature, hot working fluid temperature, temperature ratio 
and irreversible property were considered as decision variables. 
In addition, the effects of the referred decision variables on each 
objective-function like as thermal efficiency of dish-Stirling system 
were carried out while different decision making approaches 
implemented. Due to obtained solutions from different decision 
making approaches, fuzzy decision making has higher maximum 
error percent than other two addressed decision making ap¬ 
proaches such as LINMAP and TOPSIS approaches. This fact could 
be illustrated by independency of the fuzzy decision making ap¬ 
proach on ideal solution and non-ideal solution. In other hand, 
based on maximum error percent, obtained routs of TOPSIS ap¬ 
proach has superiority than other approaches like as LINMAP. Opti¬ 
mal results which are closer to ideal result are chosen. Stirling 
engine thermal efficiency is between 40% and 45% and dish-Stirling 
thermal efficiency is between 35% and 40% in turn, the results of 
the presented thermo-economic model for dish-Stirling imple¬ 
menting NSGA-II algorithm are acceptable. It should be noted that, 
Starting with the scheme for computation of the Performances 
(Efficiency and Power) for Solar Stirling Engines developed on 
Thermodynamics with Finite Speed [6-12] future papers which 
will utilize the methods described in the present communication 
will obtain results closer to the reality. The authors intend to do 
such a progress in the near future. 
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